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ABSTRACT 

We discuss the properties of magnetically dominated jet-like outflows from stellar 
mass black holes surrounded by debris tori resulting from neutron star disruption. 
These jets may have narrow cores (along the rotation axis) which are almost free 
of baryons and attain very high bulk Lorentz factors <> 10 6 . The jets give rise to a 
characteristic MeV to TeV emission as well as to relativistic shocks producing the usual 
MeV bursts. Because the outflow is highly directional the properties of the observed 
gamma-rays would depend on the viewing angle relative to the rotation axis. Even 
for the most intense bursts, which under the assumption of isotropic emission and 
substantial redshifts would be inferred to emit 10 52 — 10 53 erg, the efficiencies required 
are only 1CT 2 - 10~ 4 . 



Subject headings: gamma-rays: bursts - black hole physics - magnetic fields 



1. Introduction 

Gamma ray bursts pose two sets of problems. The first is to account for the required 
large and sudden energy release. The prime candidate here is the formation of a compact 
object or the merger of a compact binary: this can trigger the requisite energy release (few 
10 51 erg s" 1 for an isotropic burst at cosmological distances), with a characteristic dynamical 
timescales as short as milliseconds. The second problem is how this energy is transformed into a 
relativistically-outflowing plasma able to emit intense gamma rays with a nonthermal spectrum. 
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The literature on the second of these problems is already extensive. There have, in particular, 
been detailed calculations on the behavior of relativistic winds and fireballs. We have ourselves, in 
earlier papers (e.g. Rees & Meszaros , 1992, 1994, Meszaros , Laguna & Rees, 1993, Meszaros , 
Rees & Papathanassiou, 1994 [MRP94], Papathanassiou & Meszaros , 1996; also Paczyhski, 1990, 
Katz, 1994, Sari, Narayan & Piran, 1996) addressed the physical processes in relativistic winds 
and fireballs. Motivating such work is the belief that compact objects can indeed generate such 
outflows. There have, however, been relatively few attempts to relate the physics of the outflow 
to a realistic model of what might energize it. (Among these are early suggestions by Paczynski, 
1991, Meszaros & Rees, 1992, Narayan, Paczynski &; Piran, 1992, Thompson, 1994, and Usov, 
1994). These have involved either the reconversion of a burst of neutrino energy into pairs and 
gamma rays, or else strong magnetic fields. Although the former model cannot be ruled out, the 
beamed neutrino annihilation luminosity being marginal in black holes with a disrupted neutron 
star torus (Jaroszyriski, 1996), initial calculations for neutron star mergers are discouraging (e.g. 
Ruffert et.al. 1997); we here focus attention on magnetic mechanisms. In the present paper we try 
to incorporate our earlier idealized models into a more realistic context, and consider some of the 
distinctive consequences of an outflow which is directional rather than isotropic. 



2. Magnetically Driven Outflows 

Magnetic fields must be exceedingly high in order to transform rotational energy quickly 
enough into Poynting flux. At the 'base' of the flow, at a radius r\ ~ 10 6 — 10 7 cm, where the 
rotation speeds may be of order c, the field strength must be at least ~ 10 14 G to transmit 10 51 
ergs in a few seconds. These fields are of course higher than those in typical pulsars. However 
, as several authors have noted, the field could be amplified by differential rotation, or even by 
dynamo action (e.g. Thompson & Duncan, 1993). If, for instance, a single fast-rotating star 
collapses or two neutron stars spiral together producing (even very transiently) a differentially 
rotating disc-like structure, it need only take a few dynamical timescales for the field to amplify 
to the requisite level (which is, incidentally, at least two orders of magnitude below the limit set 
by the virial theorem). 

The most severe constraint on any acceptable model for gamma-ray bursts is that the baryon 
density in the outflow must be low enough to permit attainment of the requisite high Lorentz 
factors: the absolute minimum, for wind- type models invoking internal dissipation, is T ~ 10 2 ; 
impulsive models depending on interaction with an external medium require F ~ 10 3 . Since the 
overall efficiency is unlikely to exceed 10 -1 , this requires any entrained baryons to acquire at least 
10 3 times their 'pro rata share' of the released energy. There are of course other astrophysical 
situations where an almost baryon-free outflow occurs - for instance the wind from the Crab 
pulsar, which may contain essentially no baryons. However, this is not too surprising because the 
internal dissipation in pulsars is far too low to generate the Eddington luminosity. On the other 
hand, in GRBs the overall luminosities are ^ 10 13 LecI- It is hardly conceivable that the fraction 
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channeled into thermal radiation is so low that radiation pressure doesn't drive a baryon outflow 
at some level. The issue is whether this level can be low enough to avoid excess 'baryon poisoning'. 

When two neutron stars coalesce, some radiation-driven outflow will be induced by tidal 
dissipation before coalescence (Meszaros &; Rees 1992). When the neutron stars have been 
disrupted, bulk differential rotation is likely to lead to more violent internal dissipation and a 
stronger radiation-driven outflow. Almost certainly, therefore, some parts of the outflow must, for 
some reason, be less accessible to the baryons. 

3. Axisymmetric Debris and Jets Around Black Holes 

One such reason might be that the bursts come from a black hole orbited by a disk annulus or 
torus (e.g. Paczyhski 1991, Levinson & Eichler, 1993). This is of course what happens when the 
central part of the rotating gaseous configuration has collapsed within its gravitational horizon; 
otherwise, there is no reason why material should avoid the centre - indeed, there is more likely to 
be a central density peak in any non-collapsed configuration supported largely by rotation (either 
a single star or a compact binary after disruption). 

Such a configuration could come about in two ways: 

(i) The spinning disk that forms when two neutron stars merge (e.g. Davies et.al. 
1994), probably exceeds the maximum permitted mass for a single neutron star; after 
viscosity had redistributed its angular momentum, it would evolve into a black hole 
(of 2-3 Mq) surrounded by a torus of mass about 0.1 M (Ruffert et.al. 1996). 

(ii) The system may result from coalescence of neutron star and black hole binary. 
If the hole mass is 5M , the neutron star would be tidally disrupted before 
being swallowed, leading to a system resembling (i), but with a characteristic 
radius larger by a factor two and a torus mass of ~ 1 instead of O.1M . 
Numerical simulations yield important insights into the formation of such configurations 
(and the relative masses of the hole and the torus); but the Lorentz factors of the outflow are 
sensitive to much smaller mass fractions than they can yet resolve. 

It is, however, a general feature of axisymmetric flows around black holes that the region 
near the axis tends to be empty. This is because the hole can swallow any material with angular 
momentum below some specific value: within a roughly paraboloidal 'vortex' region around the 
symmetry axis (Fishbone & Moncrief, 1976), infall or outflow are the only options. Loops of 
magnetic field anchored in the torus can enter this region, owing to 'buoyancy' effects operating 
against the effective gravity due to centrifugal effects at the vortex walls, just as they can rise 
from a flat gravitating disk. These would flow out along the axis. There can, in addition, be an 
ordered poloidal field threading the hole, associated with a current ring in the torus. This ordered 
field (which would need to be the outcome of dynamo amplification rather than just differential 
shearing) can extract energy via the Blandford-Znajek (1977) effect. In the latter the role of the 
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torus is mainly to anchor the field: the power comes from the hole itself, whose spin energy can 
amount to ~ 10 53 erg. Irrespective of the detailed field structure, there is good reason to expect 
any magnetically-driven outflow to be less loaded with baryons along the rotation axis than in 
other directions. Field lines that thread the hole may be completely free of baryons. 

4. Baryon-Free Outflows 

As a preliminary, we consider a magnetically-driven outflow in which baryonic contamination 
can be neglected. In the context of many models this may seem an unphysical limiting case: the 
difficulty is generally to understand how the baryonic contamination stays below the requisite 
threshold. But the comments in the previous section suggest that it is not impossible for the part 
of the flow that emanates from near a black hole and is channelled along directions aligned with 
the rotation axis. 

There have been earlier discussions (dating back to Phinney, 1982) of relativistic MHD flows 
from black holes in AGN. In our present context, the values of L/M are larger by at least ten orders 
of magnitude. This means that the effects of radiation pressure and drag are potentially much 
stronger relative to gravity; also, pair production due to photon-photon encounters is vastly more 
important. For an outflow of magnetic luminosity L and e ± 7 luminosity L w & L channeled into 
jets of opening angle 6 at a lower radius r\ = 10 6 r6 cm the initial bulk Lorentz factor is Ti ~ L/L w , 
and the comoving magnetic field, temperature and pair density are B[ ~ 2.5 x 10 14 Ll ) { 2 r^ 1 rQ 1 9~ 1 
G, T( ~ 2.5 x 10 10 L^ 4 r- 3/4 r 6 - 1/2 ^- 1 /2 K and tij~4x ltf 2 L%% 9, % 3 /V 3 / 2 cm" 3 (primed 
quantities are comoving). Unless >> 1 the jet will be loaded with pairs, very opaque and in 
local thermal equilibrium. It behaves as a relativistic gas which is "frozen in" to the magnetic 
field, and expands with T' oc r~ l . The lab-frame transverse field B oc r _1 , and the comoving 
B' ~ B/T. The comoving energy density (predominantly magnetic) is e' oc r~ 2 T~ 2 , and the the 
pair density is n' oc T' 3 oc r~ 3 so T oc n'/e' oc r, or T ~ Ti(r/ri). 

When the comoving temperature approaches m e c 2 /k the pairs start to annihilate and 
their density drops exponentially, but as long as the scattering depth t' t > 1 the annihilation 
photons remain trapped and continue to provide inertia, so T' oc r -1 , T oc r persists until 
T' a ~ 0.04m e c 2 ~ 17 keV at a radius r a , where r^(r a ) ~ 1. Between r\ and r a this leads to 
(rJrO ~ (T>/T>) ~ W 2 lI{ 4 T; 3/ % 1/2 9^/ 2 , with T a ~ r,(r B /r,) ^ WL^Th 1 ^^! 2 . At 
r a the adiabatic density n' aad ~ n;(r a /n) 3 ~ 4 x 10 26 cm -3 is mostly photons, while the pair 
density from Saha's law is n' a ~ 1.5 x 10 18 F;r f T 1 , the photon to pair ratio being ~ 10 8 . The 
annihilation photons streaming out from r a appear, in the source frame, as photons of energy 
around T a ?>kT' a ~ bLl^r 1 / 4 ^^" 1 / 2 MeV. 

Beyond r a the lab-frame pressure B 2 oc r -2 but the inertia is drastically reduced, being 
provided only by the surviving pairs, which are much fewer than the free-streaming photons. In 
the absence of any restraining force, T oc n'/e' a n' / B' 2 oc B 2 /n', and the gas would accelerate 
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much faster than the previous T oc r. (The pair density is still above that needed to carry the 
currents associated with the field, analogous to the Goldreich-Julian density, so MHD remains 
valid). However, the Compton drag time remains very short, since even after t' t < 1 when most 
photons are free-streaming, the pairs experience multiple scatterings with a small fraction of the 
(much more numerous) photons for some distance beyond r a . 

One can define an "isotropic" frame moving at Tj oc r, in which the annihilation photons 
are isotropic. In the absence of the magnetic pressure, the drag would cause the electrons to 
continue to move with the radiation at Tj oc r. The magnetic pressure, however, acting against a 
much reduced inertia, will tend to accelerate the electrons faster than this, and as soon as V ^ Tj, 
aberration causes the photons to be blueshifted and incident on the jet from the forward direction, 
so the drag acts now as a brake. In the isotropic frame the jet electron energy is 7 = T/Ti and its 
drag timescale is that needed for it to encounter a number of photons whose cumulative energy 
after scattering equals the energy per electron, td r ,i ~ m e c 2 / '(u p h,i<Tcrf) = (m e c 2 A7rr 2 Tf /Lctt^). In 
the lab frame this is Tj times longer, and the ratio of the drag time to the expansion time r/c 
must equal the ratio of the kinetic flux to the Poynting flux, n'jin e c 2 r 2 /[(Bf /47r)(r;/r) 2 ], where a 
is the scattering cross section and n'j is the comoving pair density in the jet. This is satisfied for 
r ~ T a at r a , and since the drag time is much shorter than the annihilation time, the pair number 
is approximately frozen, and T cx r 5 / 2 for r > r a . The upscattered photons will, in the observer 
frame, appear as a power law extension of the annihilation spectrum, with photon number index 
-2, extending from ~ 0. l2F a m e c 2 to ^ I\,m e c 2 . 

The acceleration T oc r 5 / 2 abates when the annihilation photons, of isotropic frame energy 
0.12m e c 2 (r a /r), are blueshifted in the jet frame to energies ^ m e c 2 . Their directions are 
randomized by scattering, and collisions at large angles above threshold lead to runaway 77 — > e ± 
(the compactness parameter still being large). This occurs when Yp ~ r 6 6 1 ' z , at 

r v ~ lO^Lj^r, 3 ^ 4 rg^ 2 ^ -1 / 2 cm. Thereafter the threshold condition implies V oc r 2 , until the 
inertial limit T in ~ 10 9 Lj{ 4 r, 1/4 rJ /2 -1 / 2 is reached at r in ~ 10 11 Lg{ 4 rf 3/4 r6(9~ :l/2 cm. Besides 
going into pairs, a reduced fraction of the Poynting energy may continue going into a scattered 
spectrum of number slope -2 up to Ti n m e c 2 . 

However, in an outflow of L ~ 10 51 L5i erg s -1 maintained for a time t w ~ few seconds, the 
relativistic jet would have to push its way out through an external medium, with consequent 
dissipation at the front end, as in double radio sources. Except for a brief initial transient (<C 1 s in 
observer frame) the shock will be far outside the characteristic radii discussed so far. The external 
shock will be ultrarelativistic, and slows down as r~ l l 2 . The jet material, moving with r > 1, 
therefore passes through a (reverse) shock, inside the contact discontinuity, which strengthens 
as r 1 / 2 as the external shock slows down. Since it is highly relativistic and the medium highly 
magnetized, this reverse shock can emit much of the overall burst energy on a time ^ t w . When 
t ~ t w the external shock has reached ~ l0 1Q L l J±no l ^t\i 2 6~ 1 / 2 cm, where n a is the external 
density, and the Lorentz factor has dropped to ~ lO 3 -^ 8 ?!,, l ^t w l ^ 4: 6~ 1 l A . This, as well as the 
expected radiation of the shocks at (and after) this stage is rather insensitive to whether the initial 



-6- 



Lorentz factor is indeed ~ 10 6 — 10 8 , or whether baryon loading has reduced it to ~ 10 3 . After 
this the flow is in the impulsive regime and produces an "external" shock GRB on an observer 
timescale rd/cT 2 ~ t w , as in, e.g. Meszaros & Rees, 1993, MRP 1994. 

5. Radiation from High T Magnetically Driven Jets 

If the jet is indeed baryon-free, and therefore has a Lorentz factor ~ 10 7 — 10 9 , an extra 
mechanism can tap the Poynting energy along its length (before it runs into external matter) 
- namely, interaction of pairs in the jet and annihilation photons along the jet with an ambient 
radiation field. In our context, this would be radiation emitted by the torus or baryon-loaded 
wind that is expected to flow outward, mainly in directions away from the rotation axis, forming 
a funnel that surrounds the jet. 

The ambient radiation causes an additional drag, which limits the terminal Lorentz factor of 
the jet below the values calculated in §4 in those regions where this is important. As a corollary, 
the Poynting flux is converted into directed high energy radiation, after pair creation in the jet 
by interaction with the annihilation radiation. We cannot generally assume that the ambient 
radiation field is uniform across the whole jet, because it may not be able to penetrate to the axis, 
but the boost in photon energy can in principle be ^ T 2 . This is similar to what is discussed in 
AGN but with more extreme Lorentz factors and radiation densities. Since the ambient radiation 
has a luminosity *5 Lecl-i an d the burst luminosity is larger by 10 12 — 10 13 , this mechanism is 
significant only when the jet Lorentz factor has the very high values <> 10 6 characteristic of these 
baryon-free outflows. 

The photons from the sides of the funnel emitted into the jet have energies 
Xf = i? 7 /m e c 2 ~ 1/20. When they pass transversely into the jet, and are scattered by 
"cold" pairs, the scattering will be in the K-N regime for any F > 20. The electron will therefore 
recoil relativistically, acquiring a Lorentz factor ~ T/20 (in the comoving frame), and will then 
cool by synchrotron emission, isotropic in the jet frame. This radiation will, in the observer 
frame, be beamed along the jet and blueshifted by T. One readily sees that the net amplification 
(energy/incident energy) is ~ T 2 . If the external photons instead interact with one of the beamed 
gamma-rays of energy E in the source frame (typically near threshold for pair production) the 
resultant pair will have a Lorentz factor of order T/(E/m e c2) in the comoving frame, and will 
again emit synchrotron, yielding the same amplification factor as before. 

The interaction of the ambient photons with the beamed annihilation photons 
in the jet leads to a complicated cascade, where the jet Lorentz factor must 
be calculated self-consistently with the drag exerted by the ambient photons. 
A schematic cascade at r ~ r p where T p ~ 10 7 , would be as follows. 
A) Interactions near r p between funnel photons of lab frame energy xj ~ 10 _1 and 
beamed annihilation jet photons of lab frame energy x\ ~ x a ~ 10 (produced by pair 
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recombination in the outflow) lead to of energy 72 ~ x\ ~ 10 (lab), or 72 ~ r p /72 ~ 10 6 
(comoving jet frame). In the comoving magnetic field B 1 ~ 6 x 10 3 G these pairs produce 
synchrotron photons of energy x' 2 ~ 10 1 (comoving), or x 2 ~ x 2 r p 2 x 10 8 (lab). 

B) Photons %2 interact with ambient photons Xf to produce pairs of energy 

73 ~ X2 ~ 2 x 10 8 (lab), or 73 ~ 73 /Tp ~ 2 x 10 1 (comoving). In the same comoving 
field these produce synchrotron photons x' 3 ~ 1CP 8 (comoving) or x 3 ~ 10 _1 (lab). 

C) Photons X3 ~ 10 _1 are below threshold for pair production with funnel photons Xf ~ 10 -1 , 
ending the cascade. The resulting photons have lab energes E% ~ xsm e c 2 ~ hOLlJ-^T^rQ l l 2 Q~ l l 2 
KeV. Of course, in a more realistic calculation the self-consistent jet Lorentz factor may vary 
across the jet, and one needs to integrate over height. However this simple example illustrates the 
types of processes involved. 



6. Discussion 

lYiom. the scenario above it follows that Poynting dominated (or magneticaly driven) outflows 
from collapsed stellar objects would lead, if viewed along the baryon-free jet core, to a 7-ray 
component resembling the GRB from the external shocks in 'impulsive fireball' models. However 
the reverse shock is here relativistic, and may play a larger role than in impulsive models. Viewed 
at larger angles, the GRB component would ressemble that from internal shocks (see below). The 
characteristic duration t w and the (possibly complicated) light curve are controlled by the details 
of the Poynting luminosity production as a function of time, e.g. the cooling or accretion time of 
the debris around a centrally formed black hole. 

Besides the 'standard' GRB emission, Poynting dominated outflows viewed along the 
jet core would be characterized by additional radiation components from the jet itself (§§4, 
5). The annihilation component peaks at E x ~ SL^r^ 2 ^ 17 V 1 / 2 MeV (or 50 MeV if 
9 ~ 0.1, T[ ~ 10), with a power law extension of photon number slope -2 going out to 
& T p m e c 2 ~ SLg^r^rg 1//2 -1 / 2 TeV, if ambient photon drag limits r to & T p (otherwise, it could 
extend to Ti n m e c 2 ~ 500 TeV). If, as argued by Illarionov & Krolik (1996) for AGN, an outer 
skin of optical depth unity protects the core of the jet from ambient photons, this annihilation 
component could have a luminosity not far below the GRB MeV emission, and the T of the jet core 
would not be limited by ambient drag, only that of the skin. However, the skin depth is hard to 
estimate without taking the drag self-consistently into account. The cascade from the interaction 
of ambient (funnel) photons with jet photons leads to another jet radiation component. The 
simplified estimate of §5 gives as a characteristic energy, with a 

power law extension. The amplification factor of the cascade is A c ^ T 2 ^ 10 14 Lg^r^r^ 1 ^ 1 , 
and since the funnel emits ~ LEdd, the luminosity could be a (possibly small) fraction of L. The 
duration of these components is ~ t w , preceding the normal MeV burst by ~ t w , and it may be 
more narrowly beamed, arising from regions with T p ~ 10 7 , as opposed to T & 10 3 for the MeV 
components, and T a ~ 10 2 for the peak annihilation component. 
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Gamma-rays will be detected for a time t w (longer if external shocks are expected). There 
may also be prolonged after effects, e.g. Meszaros & Rees, 1997, with X-ray and optical fluxes 
decreasing as a power law in time. In a rotating black hole - torus system with a super-Eddington 
outflow the baryon contamination will be minimal along the rotation axis and will increase at 
larger angles to it. A narrow, largely baryon-free jet core such as discussed in §§4 and 5, would 
be surrounded by a debris torus producing a baryon-loaded, slower outflow acting as a funnel, 
which injects X-ray seed photons into the jet. This slower outflow would carry a fraction of the 
GRB luminosity L in kinetic energy form. For a slow outflow T ~ 10, this kinetic energy can be 
reconverted into nonthermal X-rays after a time t x ~ day if it is shock-heated at a radius r ~ 10 15 
cm, either by Alfven wave heating, or by interaction with a pre-ejected subrelativistic shell of 
matter of ~ 1O _3 M as has been detected in SN 1987a. This could lead to the substantial X-ray 
afterglow detected ~ days later in GRB 970228 (Costa, et al., 1997). 

The observed time structure and luminosity of the gamma-ray burst would depend on the 
angle between the rotation axis and the line of sight. Viewed obliquely, the outflow has T ^ 10 
and only an X-ray transient with some accompanying optical emission would be seen. At smaller 
angles to the jet axis, outflows with V & 10 2 would also be seen (which might originate at r > 10 10 
cm by entrainment of slower baryonic matter by the faster jet, or they might already originate 
closer in), and the predominant radiation would arise from internal shocks (Rees Sz Meszaros , 
1994, Papathanassiou & Meszaros , 1996), which can have complex, multiple-peaked light curves. 
Closer to the rotation axis, outflows with T ~ 10 3 may dominate, either at large radii or already 
lower down, and radiation from external (deceleration) shocks would be prominent (e.g. MRP94). 
Nearest to the rotation axis, if there is indeed an almost baryon-free core to the jet, the annihilation 
radiation power law component, the relativistic reverse shock radiation and the cascade process 
(§§4,5) yield significant extra contributions, which arrive ahead of the external shock burst. The 
luminosity function for the burst population may in part (or even entirely) be determined by the 
beam angle distribution, jet edge effects and the observer angle relative to the jet axis. A time 
variability would be expected from intermittency in the Poynting flux extraction, or wobbling of 
the inner torus or black hole rotation axis, on timescales down to t v ~ ri/c ~ 10 -3 — 10~ 4 s. 

If the Poynting flux is derived from the accretion energy of a NS-NS remnant torus of 0.1 
Mq, and the gamma-rays are concentrated within an angle 6 ~ 10°, the efficiency of conversion 
of rest mass into magnetic energy need only be 10 -4 to generate a burst which (if isotropic) 
would be inferred to have an energy of 10 51 ergs, and only 10~ 2 to simulate a 10 53 erg isotropic 
burst. For a BH-NS merger the torus may be ~ 1M , and the corresponding magnetic efficiencies 
required are 10~ 5 and 10~ 3 respectively. Thus, even bursts whose inferred isotropic luminosites 
are 10 52 — 10 53 erg, either due to high redshifts or exceptional observed fluxes, require only very 
modest efficiencies. Further reductions in the required efficiencies are possible if the Poynting flux 
is due to the Blandford-Znajek mechanism, which can extract ^ 10 _1 of a near-maximally rotating 
Kerr black hole rest mass (as would be expected from a NS-NS merger), leading to equivalent 
isotropic energies ~ 10 53 (47r/27r# 2 ) 3> 10 53 erg. Even without the latter, somewhat speculative, 
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possibility, it is clear that jet-like Poynting flows from tori around black holes can produce even the 
most intense bursts with comfortably low efficiencies. The detailed burst properties - particularly 
the rapid variability and bursts seen along the jet core, the blueshifted annihilation and high 
energy cascade photons- can help to pin down the parameters of the model, while the afterglow 
can provide information on the dynamics and mass flux in directions away from the jet axis. 

We thank NASA NAG5-2857, NATO CRG-931446 and the Royal Society for support, the 
Institute for Advanced Studies, Princeton, for its hospitality, and Ira Wasserman for useful 
comments. 
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